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Abstract 

We review recent developments on axion cosmology. Topics include : axion cold 
dark matter, axions from topological defects, axion isocurvature perturbation and 
its non-Gaussianity and axino/saxion cosmology in supersymmetric axion model. 
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1 Introduction 

The standard model (SM) in particle physics is well established after the discovery of 
the Higgs-boson like particle at LHC [H [2]. While it goes without saying that SM is a 
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successful theory, it suffers from the strong CP problem [3|, H] . The following term in the 
Lagrangian allowed by the gauge symmetry, 

violates CP and contributes to the neutron electric dipole moment (NEDM). Here Qs 
is the QCD gauge coupling constant, G^^, is the gluon field strength, G^^, is its dual, 
and 6' is a constant parameter. Recent experimental bound on the NEDM reads \dn\ < 
2.9 X 10~^^ecm (90% CL) [5]. This leads to the constraint on the 6 parameter as 6' < 
0.7 X 10"^^ [21 II]. There seems to be no reason in the SM why 6 must be so small : this 
is the strong CP problem. 

Peccei and Quinn [6l [7] proposed a beautiful solution to the strong CP problem. 
They introduced anomalous global U(l) symmetry, which we denote by U(1)pq called PQ 
symmetry, which is spontaneously broken. Then the 6 term is replaced by a dynamical 
field, which automatically goes to zero by minimizing the potential. It was soon realized 
that such a solution to the strong CP problem leads to a light pseudo scalar particle : ax- 
ion [HI [9] . The axion is a pseudo-Nambu Goldstone boson in association with spontaneous 
breakdown of the PQ symmetry. It has a coupling as 

where a denotes the axion field and Fa is the scale of PQ symmetry breaking. Then 
the 6 parameter is effectively replaced with 6 + a/ Fa- It has a CP-conserving potential 
minimum at 6 + a/Fa = 0, hence the CP angle is dynamically tuned to be zero without 
fine-tuning. It is a very attractive idea for solving the strong CP problem. 

Once we believe the PQ solution to the strong CP problem, the axion may play an 
important role in particle phenomenology and cosmology. In this article we review the 
axion cosmology, in particular focusing on recent developments in the last few years. We 
do not aim to explain underlying physics of the axion, for which we refer to excellent 
reviews [31 H]. In Section [2l the PQ and axion models, and experimental/observational 
constraints are briefly summarized. In Section |3l we discuss the axion cosmology. In 
particular, recent calculations on the cold and hot axion abundances, axions emitted 
from topological defects, and axion isocurvature fluctuation and its non-Gaussianity are 
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summarized. In Section HI we focus on supersymmetric (SUSY) axion models and their 
cosmological effects. There are some recent developments on the evaluation of the saxion 
and axino abundances. In Section [5], we mention some related topics which are not covered 
in the main text. 

2 Peccei-Quinn mechanism and axion 
2.1 Models of invisible 2Lxion 

Early models of axion [8l [9] , where the axion was associated with the weak scale Higgs 
boson, were soon ruled out experimentally. Currently the most axion models make the 
axion invisible by assuming very high PQ scale. There are two known class of invisible 
axion models : KSVZ model (or also called hadronic axion model) [lOl [11] and DFSZ 
model [T2l[T3]. 

In the KSVZ model [lOl [11], heavy quark pair, Q and Q are introduced which are 
fundamental and anti-fundamental representations of SU(3)c and couple to the PQ scalar 
(j) as 

C = k(f)QQ. (3) 

Here U(1)pq charges are assigned as0(+l), Q{—l/2) and Q{—l/2). These quarks become 
heavy after gets a vacuum expectation value (VEV) of r], which is related to the PQ 
scale Fa defined in (|2]), through the relation Fa = rj/N-o^^. Here a model-dependent integer 
iVow is called the domain wall number (See Section [3]). All the SM fields are assumed to 
be singlets under U(1)pq. Clearly, this global U(1)pq has an anomaly under the QCD. 
Therefore, the axion obtains a coupling as Equation [2] and the theta angle is dynamically 
tuned to be zero by the PQ mechanism. For a minimal case where only one pair of heavy 
quarks is introduced, we have A^dw = 1- There is no domain wall problem in this case. 
(See Section [3]). 

In the DFSZ model [121 [IS], the PQ field couples to the SM Higgs. In this model 
two Higgs doublets are required : Hi and H2- We assume that Hi transforms as the SM 
Higgs and H2 as its conjugation under the SM gauge groups. Moreover, PQ charges are 
assigned so that the combination H1H2 has —1. Then we can write down the interaction 
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term among them in the potential of PQ and Higgs sector, 

- £ = |0|2(ci|i7i|2 + csli/sH + cs\H^\''\H2\^ + c^lH^H^l'' + {^l(|)H^H2 + h.c), (4) 

where ci, . . . , C4 are numerical constants and /i is a dimensional parameter. Let us suppose 
that the VEVs of Hi and H2 give up- and down-type quark masses. Then up-type SM 
quarks necessarily have PQ charge and the U(1)pq becomes anomalous under the QCD. 
Hence the axion coupling like Equation [2] appears, solving the strong CP problem. In 
this case, the domain wall number is calculated as A^dw = 3 reflecting three family of SM 
quarkslll Hence it may suffer from the domain wall problem if the PQ symmetry is broken 
after inflation. A distinct feature of the DFSZ model is that SM fermions, including 
leptons, have tree-level coupling to the axion through the Higgs-axion mixing. 

In the so-called variant axion model [TH |15], it is assumed that Hi only couples to 
one family of up-type quarks (say, top quark). All other quarks have zero PQ charges 
and obtain masses from H2. In this case, we have A'^dw = 1 and there is no domain wall 
problem. See [16] for implications of Higgs sector in this model at collider experiments. 



2.2 Astrophysical and experimental constraints 

Since the axion interaction is very weak and much lighter than typical temperature of 
stars, axions are emitted from stars. Hence the PQ scale is bounded below so as not 
to change evolutions of stars significantly [T71 118] . The most stringent constraint comes 
from the observation of SN1987A. In order for the axion emission not to shorten the 
burst duration, the PQ scale is bounded as ^ 4 x 10^ GeV. This bound relies on the 
axion-hadron interaction, which always exists in the axion model solving the strong CP 
problemf" 

lower bound on the PQ scale as Fa ^ 10^ GeVlf 

There are activities on experimental searches for the axion. 



The observation of horizontal branch (HB) stars in globular clusters also set 



^ This depends on PQ charge assignments on the Higgs field. For example, if we take the PQ charge 
of H1H2 to be —2, the allowed term is (jpHiH2 and we obtain iVow = 6. 

^ For much smaller Fa the axion becomes optically thick and the axion emission rate is suppressed. 
The bound from the burst duration of SN1987A disappears at Fa ^ 10^ GeV. On the other hand, for even 
smaller Fa, the detection rate of the thermally emitted axion from SN1987A at the Kamiokande would 
become large. As a result, the SN1987A does not pose a constraint for 10^ GeV ^ Fa^ 10^ GeV. 

^ The bound from HB stars come from axion-photon interaction, which is model-dependent. If the 
coefficient of the axion-photon-photon interaction is somehow chosen to be much smaller than unity, the 
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Axion helioscope Axion helioscopes try to detect axions with energy of order keV 
emitted from the center of the Sun through the axion-photon conversion process under 
the magnetic field [21]. Currently the CAST experiment [22l |23] put the most stringent 
bound on the strength of axion-photon-photon coupling in this way. In particular, the 
resonant conversion expected from the plasma frequency of photon induced by the buffer 
gas can improve the sensitivity and they begin to cover the parameter region predicted 
in the QCD axion for Fa ~ 10^~® GeV. The Tokyo axion helioscope also reaches the QCD 
axion prediction for small range of the axion mass [24]. Future project, called lAXO [25] . 
may reveal the QCD axion for ~ 10^ GeV. 

Axion haloscope Axion haloscopes try to detect dark matter (DM) axions in the 
Galaxy by using the microwave cavity [211 126]. Under the magnetic field, the axion DM 
may produce radio wave with its frequency corresponding to the axion mass, and it is 
amplified if the size of the cavity matches with the Compton wave length of the axion. 
Note that this technique crucially relies on the assumption that the observed DM consists 
of cold axion. The ADMX experiment already begins to exclude axion DM for a limited 
range of the axion mass around F„ ~ 10^^ GeV [27]. It is expected to cover the wide range 
of parameters consistent with axion DM. 

Laser searches Axion mixes with the photon in the external magnetic field due to the 
axion-photon-photon interaction term. The "light shining through a wall" experiments 
utilize the laser, a part of whose light would pass through a wall as axions under the 
magnetic field. The ALPS collaboration placed the most strict bound by this idea |28j . 
although currently it does not reach the CAST sensitivity and the bound from HB star 
cooling. Ideas to significantly improve the sensitivity of laser search is proposed [29] . 

bound from HB star cooling can be avoided. This, combined with the argument from SN1987A, leads 
to the so-called "hadronic axion window" at Fa ~ 10^ GeV {ma ~ 1 eV) where constraints from stars 
are absent [El [20] . (In the DFSZ model, the axion-electron interaction induces rapid star cooling for 
this value of F^, hence there is no such an window.) However, recent cosmological data disfavor such a 
scenario, because thermally produced axions contribute as hot dark matter component of the Universe. 
(See Section EXll). 



6 



Long range forces Since axion is very light, it can mediate macroscopic long range 
forces [30]. Due to the CP-odd nature of the axion, the coherent interaction (or monopole- 
monopole interaction) between macroscopic objects are suppressed. Instead, searches 
for monopole-dipole interactions provide limits on the axion-nucleon and axion-electron 
coupling [311 [32], although it is still far from the prediction of QCD axionQ 

To summarize, the best limit on the PQ scale and the axion mass comes from astro- 
physical arguments. Future axion helioscope experiment may reach a realistic parameters 
predicted in the QCD axion model. The cavity experiment will also be sensitive to the 
axion DM : if axion is the dominant component of DM, the ADMX will detect its signa- 
tures for realistic values of Fa. Some novel ideas are also proposed for detecting the axion 
DM for F„> lO^s GeV [31]. 

3 Axion cosmology 
3.1 Evolution of PQ scalar 

The Peccei-Quinn (PQ) scalar field has the following lagrangian: 

£ = ^|9^0p-Kfr(0,T), (5) 
where ^(0, T) is the effective potential at temperature T and given by 

Kff = ^(|0p-r/T + ^T^l0p. (6) 

The above lagrangian is invariant under the global t/(l)pQ transformation, </> — </)e*" with 
a constant. At high temperature T > Tc = V^r], the potential has the minimum at = 
and the vacuum has the U (1)pq symmetry (Figure[I]). However, as the cosmic temperature 
decreases, the vacuum with = becomes unstable and the PQ scalar obtains vacuum 
expectation value |0| = r] (Figured]). Thus, f/(l)pQ symmetry is spontaneously brokem 
at T < Tc. The axion a is a Nambu-Goldstone boson associated with this spontaneous 
symmetry breaking and corresponds to the phase direction of the PQ scalar as 

= \(f)\e'^- = |0|e*'^/''. (7) 
Since the phase direction is fiat in the potential E] the axion is massless at this point. 



^ It is argued that searches for monopole-nionopole force combined with star coohng constraints give 
stronger constraint [55] . 
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Figure 1: Potential of the PQ scalar. 

When the [/(1)pq symmetry is spontaneously broken, one-dimensional topological 
defects called axionic strings are formed. After formation the axionic string networks 
evolve by emitting axions and follow the scaling solution. The emission of axions from 
the strings and its contribution to the present density is described in Section 13.3.11 

When the cosmic temperature further decreases and becomes comparable to the QCD 
scale A (~ 0.1 GeV), the axion obtains its mass through the QCD non-perturbative effect. 
Then, the axion potential is written as 



= T?r- 1 1 - cos — ) , (8) 



where is the mass of the axion. Here A^dw is called the domain wall number which is 
a model-dependent integer related to the color anomaly. The axion mass depends on 
the temperature and it was first calculated in References |36]. More recently, Wantz and 
Shellard [37] obtained the temperature dependence of nia using the interacting instanton 
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liquid model [3H]. They gave a simple approximation in the form of the power-law as 

|4.05xl0-|(f)-"' r>0.26A. 
' \ 3.82 X 10-2i; T < 0,26A, ^ ' 

where A ^ 400MeV and Fa is the axion decay constant given by Fa = tj/Ndw- The 
axion potential [S] explicitly breaks f/(l)pQ symmetry into its subgroup Z^^^ and has 
degenerated minima (vacua) at a = 0, 27r/A^DW) ■ " " 27r(A''Dw ~ 1)/^dw- When the Hubble 
parameter H becomes comparable to the axion mass rria, the axion starts to roll down to 
one of the minima. Since the axion field settles into different minima in different places of 
the Universe, domain walls are formed between different vacua. These domain walls attach 
the axionic strings. For A'^dw = 1 the domain walls are disk-like objects whose boundaries 
are axionic strings and collapse by their tension. On the other hand, the string-domain 
wall networks are very complicated for A^dw > 2 and dominate the Universe soon after 
their formation, which causes a serious cosmological difficulty called domain wall problem. 
More details about the axion domain walls are found in Sections 13.3.21 and 13.3.31 

3.2 Abundance of the axion 
3.2.1 Cold axion 

The axion field starts to oscillate when the cosmic expansion rate H becomes comparable 
to the axion mass m^. Using Equation [9l the condition ma{Ti) = 3H[Ti) leads to 

Ti = 0.98 GeV ( " ) ( ). (10) 

V10i2GeV/ V400MeVy ^ ^ 

The oscillation of the coherent field a is described by the equation of motion, 

R 

'd + 3-d + m'^{T)a = 0, (11) 
R 

where R is the scale factor. From Equation [TT] we get 

Pa = ~3^a^ + TfiaTnaa^, (12) 

where Pa{= a^/2 -|- m^a^/2) is the axion energy density. By averaging the above equation 
over an oscillation ((d) = pa and {m?^a?) = pa), it is found that paR? /rria is invariant under 
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adiabatic condition ma ^ H. Thus, we obtain the present axion number to entropy ratio 
as 

yicoM) = !hlfi=p( , (13) 

So V "5 /t=Ti 

where s is the entropy density and Sq is its present value. Here /3 is the correction 
factor taking into account that the adiabatic condition {nia ^ H) is not satisfied at the 
beginning of the oscillation. The correction factor was calculated by |39] which gives 
/3 = 1.85. Thus, the present axion density is given by 



Qah^ = 0.18 ( J" V''' ( ) , (14) 

Hl0i2GeVy V400MeVy' ^ ^ 

where h is the present Hubble parameter in units of lOOkm/s/Mpc. Here 6i = ai/t] is the 

initial angle at onset of oscillation. When the PQ symmetry is spontaneously broken after 

infiation, 6i is random in space and hence we should replace 6f by its spatial average, i.e. 

(Of) = 71^/3 X Canh, whcre Canh(— 2) is the anharmonic correction |10l E]. On the other 

hand, if PQ symmetry is broken before or during infiation, 6i takes the same value in the 

whole observable Universe. Then, 6i is considered as a free parameter. 

The density of the coherent axion oscillation cannot exceed the present DM density of 

the Universe determined from the observations of cosmic microwave background (CMB), 

ficDM^^ = 0.11. This gives the following upper bound on the axion decay constant: 

Fa < 1.4 X 10^^ GeV, (15) 

when the PQ symmetry is broken after infiation. For the case of PQ symmetry breaking 
before or during infiation, see Section 13.31^ 



3.2.2 Hot axion 

Axions are also produced in high-temperature plasma [l3l HI] . The abundance of such 
hot axions in the KSVZ model, in terms of the number-to-entropy ratio Ya = ria/s, was 
estimated recently in 



^ Notice that Eg uat ion [HI assumes no late-time entropy production after the QCD phase transition. If 
there is a late-time entropy production by decaying particles, the abundance is reduced and upper bound 
on the PQ scale is relaxed [^ . 
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where Tr denotes the reheating temperature after inflation. Notice that the above ex- 
pression apphes for Tr < with 



Otherwise, axions are thermahzed at T > and decouple at T ^ To. Then the relic 
abundance is given by = 0.28/^,(rD) 2.6 x IQ-^ for ^*(rD) = 106.75). 

The above calculations were carried out in quark-gluon plasma. For small -Fa(~ 10^ GeV), 
the axion decoupling may occur after the QCD phase transition. The decoupling tempera- 
ture and the resultant axion abundance in such a case were estimated in [191 SS] by taking 
account of the axion-pion interaction. Since the axion mass is given by ma^0.6eV for 
Fa'^ 10^ GeV, it contributes to the hot dark matter and such a contribution is restricted 
from cosmological observations. According to recent results |16], the constraint reads 
rria < 0.91 eV assuming massless neutrinos and rria < 0.72 eV after marginalizing over the 
neutrino mass. This closes a hadronic axion window (see Section 1?^ . 

3.3 PQ symmetry breaking after inflation 

Cosmological consequences of axions are different depending on whether the PQ symmetry 
is broken after inflation or not. When the symmetry breaking takes place after inflation, 
topological defects like strings and domain walls are formed in the course of evolution of 
the PQ scalar as shown in Section ISTTl On the other hand, if the PQ symmetry is broken 
before or during inflation, the produced strings are diluted away and the field value of 
the axion is the same in the whole observable Universe. Thus, the axion settles into the 
same minimum when it acquire the mass at the QCD scale, and hence domain walls are 
not formed. So no domain wall problem exists. However, in this case, the axion (which 
already exists during inflation) obtains large fluctuations and produces isocurvature den- 
sity perturbations which are stringently constrained by the CMB observations. In this 
section, we first consider the case where the PQ symmetry breaking occurs after inflation 
and see cosmological consequences the axionic strings and domain walls. 




(17) 
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3.3.1 Axionic strings 

When the U{1)pq symmetry is spontaneously broken, one- dimensional topological defects 
called axionic strings are produced. Since the U{1)pq is a global symmetry, the axionic 
string is a global one. Unlike local strings which eventually lose their energy by emitting 
gravitational waves, the emission of the Nambu-Goldstone bosons, i.e. axions is a domi- 
nant energy loss process for axionic strings. After formation the axionic string networks 
evolve by emitting axions and follow a scaling solution in which the energy density of the 
string networks pstr is written as 

Pstr = (18) 

where ^ is the length parameter which is constant in the scaling regime and fi is the string 
tension given by 

/z = vrF^n(^), (19) 

where 5s = 1/ VXrj is the width of the strings. ^ represents average number of infinite 
strings in a volume t'^(~ horizon volume) and is determined by numerical simulations |171 
HHl HH ISnj. Figure |2] shows the recent simulation [50] and the evolution of the length 
parameter ^ is shown in Figure El For the axionic string networks, ^ ^ 0.7 — 1.0. 

The axions emitted from the axionic string networks give a significant contribution to 
the present matter density of the Universe. This was first pointed out by Davis [51] [see 
also [SUES]]- The axion density due to the axionic strings depends on the energy spectrum 
of the emitted axions. Davis and co-workers claimed that the energy spectrum has a sharp 
peak at the horizon scale [511 [521 [53] . On the other hand, Sikivie and co-workers [5 ^ [55 1 156] 
insisted that the spectrum is proportional to 1/A; [k: axion momentum). Because the 
present axion density is given by pa = ma{0)na, it is crucial how many axions are produced 
from the strings. Since smaller number of axions are emitted for the 1/k spectrum, the 
present density of the axions is less important than that for Davis's spectrum. This 
controversy was solved by field theoretical lattice simulations [HI [50] which showed that 
the spectrum is sharply peaked around the horizon scale and exponentially suppressed at 
higher momenta as seen in Figure [H Using this spectrum, the mean reciprocal comoving 
momentum of the emitted axions (l/k) which is important in calculating the axion number 
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Figure 2: Evolution of the axionic string networks from the simulation [50]. The panels 
show the time slices at t = At^, 9tc, IQtc and 25tc where is the cosmic time corresponding 
to the critical temperature T^. The spatial scale shows a comoving length in unit of the 

horizon siz6 8bt tend — 25tc. 



density is estimated as [50] 

{l/k)R{t) - 0.23 (20) 

where (27r)/t is the momentum corresponding to the horizon scale. Thus, the present 
axion density due to the axionic strings is estimated as 

"""''' = 2°<I54v)"'(4o4v)- 

Notice that the string contribution is larger than the coherent oscillation given in Equa- 
tion [H 
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Figure 3: Time evolution of the length parameter C, from the simulation [50] . 



3.3.2 Axionic domain walls with A'^dw = 1 

Around the QCD scale, the axion acquires mass through QCD instanton effect and the 
axion potential is given by Equation [Hi Since the potential has A^w discrete minima and 
the axion settles down to one of the minima, the Universe is divided into many domains 
with different vacua and domain walls are formed separating those domains. Around an 
axionic string, the phase of the PQ scalar (6*^ = a/t]) rotates 27r, so Aqw domain walls 
attach each string. The cosmological evolutions of the string-domain wall networks are 
quite different between Adw = 1 and Adw > 2. 

First, let us consider the case of Adw = 1- The domain walls with Adw = 1 are disk- 
like objects bounded by strings and collapse by their tension. Thus, the string-domain 
wall networks are unstable for Adw = 1 and they decay soon after formation. This was 
confirmed in 2-dimensional [SZ] and 3-dimensional [5B] numerical simulations. In Figure |S] 
the evolution of the domain walls is shown and it is seen that the string-domain wall 
networks decay and disappear in the Universe. So in this case there is no domain wall 
problem. 
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Figure 4: Differential energy spectrum of radiated axions between ti = 12.25tc and t2 = 
25tc- Green (red) bars correspond to statistical errors alone (statistical and systematic 
errors). The comoving wavenumber k is in units of {2tc)~^- Note that the horizon scale 
corresponds to /c ~ 3. The scale in y-axis is arbitrary. 



However, a large number of axions are produced in the decay of the string-domain 
wall networks and they give a significant contribution to the present matter density of 
the Universe. Similarly to axions from the axionic strings, there has been controversy 
about the spectrum of the axions emitted from domain walls. In [59] it was found that 
the spectrum has a peak around which is the width of the domain walls, so axions 
from domain walls are mildly relativistic. On the other hand, the authors in [60| |6T] 
claimed that the mean axion energy is larger than that in [59] by a factor 20. This 
conclusion comes from the reasoning that the energy of the domain walls is converted 
into strings which emits axion with spectrum proportional to following [5^ |55| [56] . 
The recent 3-dimensional lattice simulation [58] settled this problem and showed that the 
axions produced in the string-wall decay are mildly relativistic with mean momentum 
{k)/R ~ Snia which is consistent with ^9]. The energy density of the axions today is 
then estimated as 
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Comparing to the contributions from the coherent oscillation (Equation [14] ) and strings 
(Equation [2T|) . it is found that the axions from the string- wall networks gives a dominate 
contribution to the present DM density. 

The total density of the cold axions is given by 




where we take ^ = 1.0 ± 0.5 in Equation [2T1 Thus we obtain the following constraint on 
the axion decay constant: 

Fa < (2.0 - 3.8) X 10^° GeV. (24) 
3.3.3 Axionic domain walls with A^dw > 2 

The wall-string networks with A^dw > 2 have complicated structures and do not decay 
contrary to the walls with A'^qw = 1- After formation, the long-lived wall-string networks 
evolve into the scaling regime as shown in the simulations [621 [57] . The domain wall 
density is written as 

Pwall = A^Y^, (25) 

where cXwaii is the wall tension given by (T„aii = 9.23maF^. Here A is the surface parameter 
which becomes constant in the scaling regime. Since the total cosmic density decreases 
as 1/t^, the domain walls soon dominate the Universe, which conflicts the standard cos- 
mology. 

To solve the domain wall problem, one can introduce an additional term in the axion 
potential which explicitly breaks the Z^^^ symmetry, 

6V = -Er]^{(j)e'^ + h.c.), (26) 

where S is a small parameter describing the size of the explicit Z^^^ breaking and S 
is the phase. This term is called a "bias" and lifts the degenerated vacua so that the 
potential has a unique minimum at a ~ 0. The energy difference between the neighboring 
vacuum produces a volume pressure and makes the domain walls accelerate toward the 
false vacuum. Thus, the false vacuum regions shrink and the domain walls annihilate 
each other. Let us estimate the decay epoch of the string-wall networks with bias. The 
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Figure 5: Evolution of the string-domain wall networks for A'dw = 1- The white lines 
correspond to the position of strings, while the blue surfaces correspond to the position 
of the center of domain walls, "tau" in each panel is the conformal time which is related 
to the cosmic time as tau = y^. 
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pressure due to the bias term is roughly given by the difference of neighboring vacua which 
is pv ~ 47rS?7^/A^DW- On the other hand, the pressure from the surface tension which 
straightens the wall up to the horizon scale is px ~ o'waii/i- The domain walls collapse 
when Pv ~ pr, so the decay time tdec is estimated as 

tdec - ^4^- (27) 



Using (Twaii ~ "^a'7^/^DW Write the above equation as 

tdec = "^^4-^7 — ' (28) 

where a ~ 18 is determined by the numerical simulation [57] . 

The string- wall networks continuously emit axions until they decay at tdec- The spec- 
trum of the emitted axions has a peak around the momentum scale determined by the 
mass of axion as shown in Figure [6] and the mean momentum of the axions is estimated 
as ea = {k)/R c:^ ^rria. This is the same as the case for Nj^^^ = 1. Then we can obtain the 
present density of axions from the string-wall networks, 

fi,..n/.^ - 0.3 AN^i!' (^) (y^i^) • (29) 

Since the long-lived walls produce more axions than the wall with A^dw = 1? the cosmo- 
logical constraint is more stringent as seen below. 

Gravitational waves are also produced from the string-wall networks. Since the char- 
acteristic length scale of the domain walls is the Hubble time ~ t, the quadruple moment 
Qij of the system is estimated as Qij ~ Pwaot^i^ ~ cfwaa-Af^. Using the quadrupole formula 
the energy of the gravitational waves is E^^ ~ G J\? al^.^^i^'^ . Thus we can write the energy 
density of the produced gravitational waves as 

where egw is the efficiency of gravitational weaves and egw ~ 5 from the simulation [63] . 
The present density of the gravitational waves is calculated as 

n^U^^l,-^^A^N,U-J^)-\^)-\ (31) 
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The amplitude of gravitational waves produced by the domain walls is too small to be 
detected in the current and future-planned experiments for allowed values of H and Fa 
(see below). 

Let us examine whether we have a consistent scenario with A^dw > 2 or not. First the 
cosmic density of the cold axions should not exceed the CDM density measured by the 
CMB observations which gives 

^a,tot = + ^a,str + ^a,wall < 0.11. (32) 

Second, the bias term shifts the minimum of the axion potential by 

(a)_ 2SiV^^F>in^ 

T] m2 + 2SiV2^F2cos5' ^ ^ 

which violates CP. From the experimental upper bound on NEDM [5], 6 should satisfy 

^ < 0.7 X (34) 

which leads to the constraints on S and Fa through Equation [33J Finally, we have the 
astrophysical constraint from SN1987A, > 4 x 10^ GeV. In Figure [7] these constraints 
are shown. If the phase of the bias term 6 is 0(1), there is no allowed region and the 
allowed region appears for 5^ 10"^. Therefore, when the PQ symmetry is broken after 
inflation, the axion models with A^dw > 2 are excluded unless the phase of the bias term 
is fine-tuned. 

3.4 PQ symmetry breaking before/during inflation 

If the PQ symmetry is broken before or during inflation, inflation expands a domain with 
some PQ phase 6a to the size larger than the present observable Universe. Thus, there are 
no topological defects and hence no axions from the axionic domain walls and strings in 
the Universe. The present axion density comes from the coherent oscillation and is given 
by Equation [H] : 

where 6a is constant in the whole observable Universe and considered as a free parameter. 
This leads to the upper bound on the axion decay constant. 

Fa < 6.6 X 10^^ 6a^-^^ GeV. (36) 
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Figure 6: The spectra of axions produced by string-wall networks for various values of 
iVDw in the numerical simulation [63]. The spectra shown in this figure are the difference 
of power spectra evaluated at two different time steps ti = A9t]~^ and t2 = 400?7~^. The 
momentum scale corresponding to the axion mass is given by /c ~ R(t2)ma ~ 2. 
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Figure 7: The various observational constraints in the parameter space of Fa and S. 
The green dashed- hne represents the parameter region satisfying ^a,tot < ^^cdm, and the 
region below this line is excluded since the relic abundance of axions exceeds the cold 
dark matter abundance observed today. The vertical dotted-line represents the bound 
from the observation of supernova 1987A. The red solid-lines represent the NEDM bound 
for 6 = 1, 10~^, and 10~^. The region above these lines is excluded since it leads to 
an experimentally unacceptable amount of CP-violation. The pink line represents the 
NEDM bound for the critical value 5 = 1.1 x 10"^. There are still allowed regions if the 
value of 6 is smaller than this critical value. In this figure, we fixed other parameters as 
iVnw = 6, ea = 1.5 and A = 2.6. 
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Notice that Fa much larger than 10^^ GeV is allowed for ^ 1. 

During inflation the axion field obtains the fluctuations whose power spectrum is given 

by 

Vsa = (37) 
Here the power spectrum Vsa is defined by 

{5a{k)5a{k')) = {2Tif5^^\k + k')'^Vsa{k), (38) 

where 5a{k) is the Fourier component of 5a{x) and (■ ■ ■) denotes the ensemble average. 
Here notice that the fluctuations of the axion field also give a contribution to the mean 
energy density. Using {5a?{x)) ~ H^^^/in'^ derived from Equations [37] and |38l Equation 
should be replaced by 



Qnh'^ = 0.18 



27rF„ 



Fa V-''f A \ 



10i2GeVy V400MeV 



) ■ (39) 



When the axion fluctuations are larger than the homogeneous value FaOa, the cosmic 
axion density is dominated by the fluctuations and the constraint on Fa and H[^f is given 
by [H] 

J/„.,<5^0X10-Gev(^)°" for f„ < ||. (40) 

On the other hand, for Fa > Hini/2n9a, Equation [36] is valid. 

The axion fluctuations turn into isocurvature density perturbations of the axions when 
the axions acquire mass at the QCD scale [611 [Ml [661 IHZl [SH [691 E] • The definition of 
the axion isocurvature perturbations Sa is 

Sa{x)=3{Ca-0, (41) 

where (a and ( are the curvature perturbations on the uniform density slices of the axion 
and the total matter, respectively. According to 6N formalism [70] [7T] the axion energy 
density Pa{x) on the uniform density slice of the total matter is given by Pa{x) = PaG^" 
where pa is the mean density of the axion. Since pa{x) oc (a^ + (5a(x))^ with = FaOa, 

ai ai 
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where al = af + {5a^). Define r as the ratio of the mean axion density pa to the total 
CDM density Pcdm, the CDM density on the uniform density slice of the total matter is 



PcT>M\x) — Pcdm 



'1 — r) + re 



Sa 



(43) 



Thus, we obtain the CDM isocurvature density perturbations S'cdm as 



CDM 



In ( 



Pcdm (a: j 



\ Pcdm / 

^ FJM , da" 
2r h r 



The energy fraction of the axion in CDM is given by 



1.6 



inf 



2ttR, 



1.19 



IQi^GeV j (400Mev) ' 



The power spectrum of the CDM isocurvature perturbations is written as 
(^cdm(^)^cdm(^')) = i'2^)'S^'\k + k')'^Vsa^,,ik), 



(44) 
(45) 
(46) 

(FA)2 + (i/inf/27r)2' ^'^^^ 
The amplitude of the isocurvature perturbations is stringently constrained by the CMB 
observations [72l[73l[711[l5l[76l[77l[78l[79]. The analysis using WMAP 7 year data [80] 
finds a = 'Pscdm(^*)/^c(^*) < O-^^ the reference scale = 0.002Mpc~^ where V( 
is the power spectrum of the curvature perturbations. This limit is translated into the 
constraint on Fa and Hi^f as 



4r' 



A;3 

(i^inf/27r)2 



9l 



inf 



V27rF„ 



inf 



\2'KFaJ VlQi^GeV 



2.38 



< 3.6 X 10-^\ 



(48) 



Here we have used Vc_{k^) = 2.43 x 10"^ 

Although the axion fluctuation 6a is Gaussian, Sa^-teim in Equation HHis not Gaussian 
and hence leads to non-Gaussianity in the isocurvature perturbations in CDM [811 [821 |83l 
EH ESI ESI EH ES]- The non-Gaussianity clearly appears in the bispectrum B{ki, k2, k^) 
of the isocurvature perturbations which is defined as [82] 



(5cdm(A?i)5cdm(4)5cdm(4)) = (27r)='55cDM(^i. ^2, k^)5^'^\ki + h + h 



BscoMiki,k2,k3) = fs 



"-1 "-2 



-^5cdm(^i)^5cdm(^2) + (2 perms) 



(49) 
(50) 
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Hubble scale during inflation Hj^, [GeV] Hubble scale during inflation Hj^f [GeV] 



Figure 8: Constraints on axion isocurvature model in the Hi^{-6a plane for Fa = lO^^GeV 
(left panel) ans Fa = 10^^GeV(right panel). Shaded regions are excluded by cosmological 
considerations. At small 6a, the constraint on Hij^f from the non-Gaussianity is slightly 
better than one from the power spectrum. When ifinf/27r > Fa the PQ symmetry restore 
during inflation. 

where fs is the non-linearity parameter which is given by l/(2r) for the isocurvature 
perturbations produced by the axion fiuctuationsE 

At present CMB is the best probe of non-Gaussianity in isocurvature perturbations 
and it was investigated by using WMAP 3 year data [89] and more recently WMAP 7 
year data [90]. The resultant 2a constraint is 

Jn'l'I = l«V5| < 140, (51) 
from which we obtain the following constraint: 



f Hi 



inf 



V27rF„ 



21 / rj X 4 . 17 . 3.56 



/ -f^inf \ f Fa 



\2TTFaJ VlOi^GeV 



< 2.5 X 10^^^ (52) 



In Figure El the cosmological constraints on axion models are shown in the ifjnf — 6a plane. 

4 Super symmetry and axion 

In this section we give a brief review on supersymmetric extension of the axion model. 
Since SUSY naturally solves the gauge hierarchy problem, we have a good motivation to 



^The parameter fs is related to f^f^ in [55] through the relation f^f^ = a'^fs- This is further 
related to /^f ^ in m through = (6/5)/^^^- 
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consider the SUSY axion model. In extending the axion model into a SUSY framework, 
there are some non-trivial features which must be addressed carefully. First, in SUSY, 
there exists a flat direction in the PQ scalar potential in addition to the massless ax- 
ion. This is due to the holomorphic property of the superpotential [91] : in the U(1)pq 
transformation — e*"^^^, the rotation parameter Uj can be complex. Therefore, by 
taking aj to be pure imaginary, the theory should be invariant under the scale transfor- 
mation, which means that there is a fiat direction in the scalar potential correspondingly 
and hence the PQ scalar is not stabilized. This argument breaks down under the SUSY 
breaking effect. Thus in order to stabilize the PQ scalar at an appropriate scale, we are 
forced to consider SUSY breaking and its effects on the structure of the scalar potential. 
Second, there appear very weakly coupled massive particles, saxion (scalar partner of 
axion) and axino (fermion partner of axion), which might have significant cosmological 
effects [921I93]. 

4.1 Stabilization of the PQ scalar 

In this subsection we show some explicit models of PQ scalar stabilization. We do not 
aim to make a complete list of stabilization mechanism, but rather show typical examples, 
which may be a good starting point for considering phenomenology of the SUSY axion 
model. 

4.1.1 Model A 

One of the simplest model for stabilizing the PQ scalar is described by the following 
superpotential : 

W^PQ = «:X(00-r/2), (53) 

where X is a singlet superfield with no PQ charge, and (0) is a PQ scalar with PQ 
charge (—1), and k is a constant taken to be real hereafter. The i?-symmetry under 
which X has a charge +2 and </>, (f) are singlet ensures this form of the superpotential. By 
taking into account the SUSY breaking effect, the (relevant portion of) scalar potential 
is given by 

V = ml\4>? + m?|0p + k'\<P4> - r/T + [\4>? + \^?) • (54) 
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where and m? represent soft SUSY breaking masses for and (f). PQ scalars are 
stabilized at (0) {(p) rj and {X) ~ 0. To be more precise, X obtains a VEV of 
{X) ~ 7713/2/^, where 7723/2 denotes the gravitino mass, by taking account of the constant 
term in the superpotential {Wq = m^i2Mp) to cancel the cosmological constant. In this 
model, the saxion corresponds to the fluctuation around the VEVs of (f) and cf) along the 
fiat direction 00 = 77^, and it obtains a mass of soft SUSY breaking scale, which might 
be of the order of TeV. The axino (a) corresponds to the light combination of and 0, 
which are fermionic components of and 0, and obtains a mass of 772^ = k(X) ^ 7713/2- 
(The other combination mixes with X and gets a mass of fi;7/.)0 

In order to solve the strong CP problem, we introduce a superpotential 

W^Ksvz = k(t>QQ, (55) 

for the SUSY version of KSVZ axion model, where Q and Q transforms as fundamental 
and anti-fundamental representations of SU(3)c.^ and both have PQ charges of —1/2. On 
the other hand, in the SUSY version of DFSZ model, we introduce 

W^DFSz = ^^HuH,, (56) 

where and denote the up- and down-type Higgs doublets. In this case, {HuHd) 
has a PQ charge of —2 and hence (some of) MSSM fields also have PQ charges. After 
obtains a VEV, equation gives the higgsino mass, so-called /i-term, as /i = \iff/Mp 
and it is of the order of TeV for A = 0(1) and Fa ~ 10^^ GeV. Thus the PQ scale may be 
related with the solution to the //-problem in MSSM [97] . 

4.1.2 Model B 

Let us consider the following superpotential, 

W^PQ = (57) 



Recently it was proposed that the PQ scalars in superpotential [53] can take a role of waterfall field in 
hybrid inflation, while X is the inflaton [MJ 1^51 [M] ■ A right amount of density perturbation is reproduced 
for Fa ~ IQis GeV. 

® To maintain the gauge coupling unification, it is often assumed that Q and Q are complete multiplets 
of SU(5) GUT gauge group. 
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where M denotes the cutoff scale and n{> 3) is an integer. The PQ scalars and have 
PQ charges of +1 and — n, respectively. The scalar potential after including the SUSY 
breaking effect is given by [98] 

V = -mM^ - ml\4>? + (a,^ + h.c.) + (1^1^ + n\4>?) , (58) 

where denotes the A-term contribution to the SUSY breaking potential. Here the soft 
mass squared of is assumed to have tachyonic form. Then is stabilized at 

r/=(0)~(m,M«-2)'^^""'\ (59) 

For example, for = 3 and M = Mp, we have r] ~ y/nisMp and it is of the order of 
lO^"*^ GeV for the soft mass ~ 1 TeV. The A-term also induces VEV of as (0) ~ 
A^rj/nis, which is smaller than rj if < irtg. The axino consists of mixtures of and 0, 
whose mass is given by ~ m^, as seen by substituting VEVs of and into Equation 
[571 Note that the axino can be much heavier than the gravitino if ^ ^,3/2 due to 
gauge-mediation effect (see below). 

4.1.3 Model C 

In the above two models, we introduced two (or possibly more) PQ scalar fields for 
giving rise to the scalar potential. PQ scalars are stabilized by the combination of SUSY 
potential and SUSY breaking potential. However, there may be multiple sources of SUSY 
breaking effects. In the gauge-mediated SUSY breaking (GMSB) model, the dominant 
contribution to the soft SUSY breaking mass often comes from the GMSB effect, while 
small but non-negligible gravity-mediation effect generally exists. In particular, for a 
gauge-singlet scalar such as the PQ scalar, the effect of GMSB on its potential is rather 
non-trivial because it arises at higher loop level [991 HOOl 1101] . 

Let us consider the KSVZ model and see how the scalar potential is generated through 
the SUSY breaking effect. For simplicity, we adopt the following messenger sector 

ly^ess = ^^X^^, (60) 

where X denote the SUSY breaking field and ^ and ^ are messengers fields. For k\(f)\ <^ 
A^mess, where Mmcss is the messenger scale, the mass splitting on the scalar and fermionic 
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components of Q,Q yield correction to the saxion mass. In the opposite case k\(j)\ ^ 
^mcss5 after integrating out the PQ quarks, the Kahler potential of the X field below the 
PQ scale is given by 

c = J d'ezxMm', (61) 

where the wave-function renormalization factor Zx logarithmically depends on |0| at the 
three-loop level. Thus we can estimate the PQ scalar potential as 

j -^m^' for A:|0|«M^ess 

d\<P\ "1 -(!i^|i^xplog^(^) for fc|0|»M„ ^^'^ 

where = {g'^/16'iT'^)\Fx/X\ is the typical soft mass induced by the GMSB effect. This 
potential has a negative slope, and hence PQ scalar is driven away from the origin. It can 
be stopped by other contributions. Since also feels the gravity-mediated SUSY breaking 
effect, there arises term 

l^fe-^)^m^/2|0r (63) 
The PQ scalar can be stabilized by the balance between Equations [62] and [63] |100] as 



V={m^§^J-^, (64) 
Stt^ m3/2 



if kr] > Mmess- If X dominantly breaks SUSY, i.e., \Fx\ = V3m^/2Mp, we need small k 
in order to obtain intermediated PQ scale. In this model, the axino is expected to obtain 
mass dominantly through radiative effect, which is evaluated as ^ {k"^ /167i'^)Aq, where 
Aq denotes the A-term contribution to the KSVZ coupling : £ = AqkcpQQ + h.c. with 
Q{Q) denoting the scalar component. 

One can construct a variant of above mentioned models. For example, one can make 
use of the scalar potential induced by the term [57] to stabilize the PQ scalar against the 
potential M ^QM M [M]. See also [Ml [M [107] for ideas to radiatively stabilize the 
PQ scalar. 

4.2 Axino cosmology 

The axinos are produced in the early Universe and they might cause cosmological disasters 
since the axino has a long lifetime. In the KSVZ model, the dominant axino production 
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process is scattering between gluons and gluinos through the axino-gluino-gluon inter- 
action. The axino abundance, in terms of its number-to-entropy ratio, was evaluated 
as [T081IT091 [no] 

^<...!5£.2.io-»»;(^)"(^). (65) 

where Tr denotes the reheating temperature after inflation. In the DFSZ model, there 
is an axino-Higgs-higgsino interaction. Since this is Yukawa coupling with a coupling 
constant of ~ y^jF^^ the axino production rate through the Higgs and higgsino scatterings 
increases at low temperature (T ~ /i). The axino abundance in the DFSZ model in such 
a process is evaluated recently as [1111 1112[ 1113] 

^ 10-5 { ^ ' { ^ y. (66) 

UoiiGeVy llTeVy ^ ^ 

This contribution does not depend on the reheating temperature Tr. 

Cosmological implication of the axino depends on whether it is the lightest SUSY 
particle (LSP) or not ^ If it is not the LSP, the axino decays into the neutralino whose 
abundance is limited from the observed DM abundance. If the neutralino annihilation 
cross section is fairly large, as in the case of higgsino or wino like neutralino, the axino 
decay can produce a right amount of neutralino for appropriate decay temperature of the 
axino [TT9lfT20] . 

If the axino is the LSP, on the other hand, it can also be produced by the late decay of 
heavier SUSY particles. Actually, as shown previously, the axino mass is comparable to 
the gravitino mass, or it can be much lighter. In this case the axino abundance depends 
on the sparticle abundance after their freeze out. Thus if the sparticle annihilation cross 
section is not so small (i.e., if it would reproduce the DM abundance if it were the LSP), 
no significant constraint is imposed. Axinos are also produced by the gravitinos, which 
decay into axinos and axions |121] . In this case the resulting axino abundance depends 
on the gravitino abundance, which may be produced thermally or non-thermally (e.g., by 
the decay of infiaton). If the axino is much lighter than the gravitino, the constraint on 
the reheating temperature can be relaxed. 



^ Evaluation of the axino mass is rather a delicate issue jll41 11151 1116j . See also |1171 1118] for recent 
calculations. 
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4.3 Saxion cosmology 

Next let us consider the cosmological effects of the saxion. Saxions are produced thermally 
in a similar way to the axino, but there also exists a contribution from the coherent 
oscillation. The evaluation of the latter contribution is quite model dependent. As shown 
in previous section, there are many models of saxion stabilization, and the scalar field 
dynamics should be calculated in each PQ stabilization model. We here give a typical 
behavior of the saxion dynamics, but one should keep in mind that the dynamics can be 
far more complicated in some concrete models, which would result in orders of magnitude 
difference in the estimate of saxion abundance. 

First, in the KSVZ model, the abundance of thermally produced saxions through the 
scattering by the gluons is given by |122] 

y<»..!il.2.:o-,;(^)"(^). (67) 

As far as we have recognized, no detailed calculations were performed for the saxion 
abundance created by the saxion-higgsino or saxion-Higgs interactions in the DFSZ model 
(denoted by Yj^^ hereafter), but we reasonably expect that it is same order as the corre- 
sponding axino abundance (Equation [66]) . 

As for the coherent oscillation, we first focus on the Model A in Sec. 14.11 In this case, 
the saxion feels the Hubble-induced mass in the early Universe and hence the saxion sits 
at the minimum determined by the Hubble mass if the Hubble-induced mass squared has 
positive coefficients. At if ~ rris, the saxion begins to oscillate around the true minimum 
with initial amplitude cxj. The abundance is thus given by 




where Tqs = {10/n'^g^y^^^/m^^dp. Note that it does not take into account various effects 
which would modify the abundance significantly. For example, thermal effects modify the 
(j) potential in the KSVZ model, which can change the saxion dynamics significantly [OH 
[951 1123] . Thermal corrections to the scalar potential are negligible in the DFSZ model. 
On the other hand, if the Hubble- mass squared has negative coefficient, the saxion is 
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Figure 9: The saxion abundance as a function of saxion mass for Fa = 
10^°, 10^^, 10^^ GeV from left to right. Thick blue lines correspond to (pa/s)^'"'^^ with 
ai = Fa and thin red lines to m„Yj^^ in the KSVZ model. Solid, dashed and dotted lines 
correspond to Tr = 10^*^, 10^ and 1 GeV respectively. 

driven away to ~ Mp during inflation, hence we expect ai ~ Mp in this casell^ 

The saxion abundance is given by sum of thermal contribution and coherent oscillation 
: p^/s = {pcr / s)^'-^'^'' (^i^^ + ^i''^)- Figure [9] shows it for the KSVZ model. The saxion 
abundance is bounded above from cosmological/astrophysical arguments depending on the 
saxion lifetime. Here we list relevant constraints depending on the saxion lifetime |126j . 

• To- ^ 10^^ sec : The axion produced by the saxion decay should not contribute to the 
effective number of neutrino species too much. Conversely, it can explain the recent 
observational claim of the existence of dark radiation \V27\ U5\ \12S\ \V29\ I13U] . 

• 1 sec ^ To- ^ 10^^ sec : Injected visible energy by the saxion decay should not alter 
the abundance of light elements synthesized by the big-bang nucleosynthesis (BBN). 

• 10^ sec ^ To- ^ 10^^ sec : Injected visible energy by the saxion decay should not distort 
the blackbody spectrum of CMB. 

• To- ^ 10^^ sec : Photons produced by the saxion decay should not exceed the observed 
diffuse extragalactic background photon spectrum. 

In such a case (ct^ ^ Mp), the axion isocurvature perturbation is suppressed, because it is propor- 
tional to S9/9 = S9i/9i ~ Hint/cTi |124j . On the other hand, if the coefficients are mildly positive, it is 
possible to create the axion isocurvature perturbation with extremely blue spectrum by considering the 
dynamics of saxion during inflation jl25j . 
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Figure 10: Cosmological constraints on the saxion abundance as a function of saxion 
mass nia- for Fa = 10^°GeV in the KSVZ model ((a) and (b)) and DFSZ model ((c) and 
(d)). In (a) and (c), the saxion decay into axions is assumed to be unsuppressed. In (b) 
and (d), it is assumed to be suppressed. See [126] for more detail. 

• To- ^ 10^^ sec : Injected photons should not ionize the neutral hydrogen too much, 
which would otherwise leads to too early epoch of reionization. 

• To- ^10^^ sec : The saxion abundance itself should not exceed the observed DM 
abundance. 

These set bound on the saxion abundance and reheating temperature for every mass 
range of the saxion. Figure [10] shows cosmological constraints on the saxion abundance 
for a wide range of the saxion mass. For more details on these constraints, see |126] and 
references therein. 

For Models B and C, the saxion is likely trapped at the origin = due to thermal 
effects in KSVZ-type models. In this case, thermal inflation |131] takes place because of 
the potential energy V{\(f)\ = 0). After thermal inflation, the saxion oscillation dominates 
the Universe and its decay causes the reheating. Cosmological implications of such a 
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scenario were discussed in |132l 11331 11341 11351 11361 1137] . It should be noticed that the 
domain wall number A^dw must be equal to one in order to avoid the axionic domain wall 
problem (see Section [373]) . 

5 Discussion 

We have reviewed recent developments in the field of axion cosmology and its SUSY 
version. In particular, according to recent simulations on axion emissions from axionic 
strings and domain walls, the upper bound on the PQ scale reads Fa^{2.0 — 3.8) x 
10^*^ GeV for A'^dw = 1 if the PQ symmetry is broken after infiation. There is no room 
for the case of A^dw ^ 2 unless the phase of explicit PQ breaking bias term is finely 
tuned. This gives a tight bound on the PQ model. If the PQ symmetry is broken 
before/during inflation, this constraint does not apply. Instead, the constraint from axion 
isocurvature perturbations requires low inflation energy scale. Recent investigations on 
non-Gaussianity in the isocurvature perturbation has shown that it gives comparable 
constraint to that from the analysis of CMB power spectrum. We have also reviewed 
SUSY axion models with some explicit examples of saxion stabilization. Recent studies 
on cosmology of saxion and axino are briefly discussed. 

Below we list some recent topics related to the axion physics which were not covered 
in the main text. 

Axion thermalization In a recent series of works |138[ I139[ 1140] it was pointed out 
that the cold axion thermalizes by the gravitational interaction and forms Bose-Einstein 
condensate, which also leads to drastic effects on the effective number of neutrinos through 
the effective photon-axion interaction. Currently there are little discussions on this claim 
and more investigations will be necessary (see also |141] ). 

SUSY axion model with 125 GeV Higgs After the discovery of the 125 GeV Higgs- 
like scalar boson [H |2], many models were proposed in SUSY framework. Some of them 
are closely related to the axion solution to the strong CP problem. In |104] . vector-like 
matter was introduced for raising the lightest Higgs mass while the PQ symmetry plays 
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an essential role to explain the appropriate mass of the vector-like matter and to solve 
the //-problem. In |142l 1143] , a singlet extension of the MSSM was considered for raising 
the lightest Higgs mass, where an appropriate value for the tadpole term in the singlet 
sector is explained by the PQ symmetry. Further LHC data may be able to confirm such 
a scenario. 

String axion and axiverse Finally, we want to make a brief comment on realization 
of the QCD axion in the string theory. The origin of global PQ symmetry is somewhat 
mysterious since even the Planck-suppressed PQ violating operators must be highly con- 
trolled. In string theory, such a shift symmetry of the axion-like field often appears after 
the compactification of extra dimensions in the zero-modes of the dilaton, NS-NS 2-form 
and R-R p-forms with p = 1,3 in type IIA and p = 0,2,4 in type IIB theories (see 
e.g. |144] ). These fields acquire masses from non-perturbative effects which break the 
shift symmetry. Since these effects exponentially depend on various parameters, axions 
with wide mass ranges in logarithmic scale are expected to exist, one of which may be 
the QCD axion: the so-called string axiverse scenario |145] . Here one should ensure that 
the mechanism of saxion/moduli stabilization does not give rise to the axion mass. This 
severely restricts the stabilization mechanism, and discussions on this topic are found in 
|146t 1147^ 1148^ 1149^ I15U] mostly for the type IIB theory. Once a successful saxion stabi- 
lization mechanism is identified, its cosmological effects can be discussed. It will be an 
interesting topic which has not been investigated in detail so far. 
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